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ABSTRACT 

We present sequential optical spectra of the afterglow of GRB 010222 obtained 
one day apart using the Low Resolution Imaging Spectrometer (LRIS) and the 
Echellette Spectrograph and Imager (ESI) on the Keck telescopes. Three low- 
ionization absorption systems are spectroscopically identified at Zi = 1.47688, 
Z2 = 1.15628, and = 0.92747. The higher resolution ESI spectrum reveals 
two distinct components in the highest redshift system at zia = 1.47590 and 
Zih = 1.47688. We interpret the zit, = 1.47688 system as an absorption feature of 
the disk of the host galaxy of GRB 010222. The best fitted power-law optical con- 
tinuum and [Zn/Cr] ratio imply low dust content or a local gray dust component 
near the burst site. In addition, we do not detect strong signatures of vibra- 
tionally excited states of H2. If the GRB took place in a superbubble or young 
stellar cluster, there are no outstanding signatures of an ionized absorber, either. 
Analysis of the spectral time dependence at low resolution shows no significant 
evidence for absorption-line variability. This lack of variability is confronted with 
time-dependent photoionization simulations designed to apply the observed flux 
from GRB 010222 to a variety of assumed atomic gas densities and cloud radii. 
The absence of time dependence in the absorption lines implies that high-density 
environments are disfavored. In particular, if the GRB environment was dust 
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free, its density was unlikely to exceed uh =10^ cm~^. If depletion of metals 
onto dust is similar to Galactic values or less than solar abundances are present, 
then TT-i/ > 2 x 10^ cm~^ is probably ruled out in the immediate vicinity of the 
burst. 

Subject headings: Gamma rays: Bursts — Cosmology: Observations — Galaxies: 
abundances, ISM, star clusters 

1. Introduction 

Years after the serendipitous discovery of gamma-ray bursts (GRBs) by the Vela space- 
craft (Klebesadel et al. 1973), the nature of the progenitors responsible for generating the 
initial explosion remains uncertain, while an elegant set of ideas about the production and 
evolution of the afterglow has developed (Meszaros & Rees 1997; Sari & Piran 1997). 

Aside from recent breakthroughs in X-ray spectroscopy of GRBs {e.g., Piro et al. 2000) 
a large fraction of the observational information about the environments of GRBs derives 
from optical spectroscopy and imaging of the host galaxies. Optical spectroscopy of the 
integrated light and calibrated emission lines has been used to derive the star-formation rate 
(SFR) for a number of the host galaxies (see for instance Bloom et al. 1998; Djorgovski et 
al. 1998). As a complement to spectroscopy, high-resolution optical imaging and astrometry 
have helped pinpoint two-dimensional locations for GRBs with respect to the host galaxy. 
The distribution of these GRB locations has yielded a number of statistical constraints in 
the progenitor scenarios {e.g., Bloom, Kulkarni, & Djorgovski 2000). 

A promising technique that has been applied with less success so far is the study of time 
dependence in absorption lines from metal ions (Perna & Locb 1998; Bottcher et al. 1999) 
or H2 vibrational levels (Draine 2000) that might be excited by the UV radiation generated 
during the evolution of the burst. The detection of absorption-line variability could provide 
important clues about the physical dimensions of the photoionized region and the density 
in the immediate environment of the GRB. Absorption-line variability can be quantified 
by measuring changes in the equivalent widths as a function of time or by identifying the 
appearance of new absorption features. Vreeswijk et al. (2001) studied the time evolution of 
the Mg II doublet in the optical spectra of GRB 990510 and GRB 990712 but failed to find 
any significant changes. In this work we present a study of the time evolution of absorption 
systems in the optical spectrum of GRB 010222, and discuss possible implications for its 
progenitor environment. The outline of the paper is as follows: §2 describes the optical 
spectroscopy, §3 describes the absorption line identification and continuum fitting. In §4 
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we detail column density determinations, arguments in favor of the z — 1.47688 redshift for 
GRB 010222, kinematics of the host galaxy, and a study of time evolution of absorption 
lines. A description of the photoionization code and results is given in §5 and §6. Finally, 
the implications of our results and conclusions are presented in §7 and §8. 

2. Observations 

2.1. LRIS Spectroscopy 

GRB 010222 was initially detected by BeppoSAX on UT 2001 Feb. 22.30799 (Piro 2001). 
An optical counterpart was reported only a few hours later by Henden (200 la, 200 lb). Our 
group began spectroscopic observations at the position of the optical counterpart on UT 
2001 February 22.66 with the dual-beam Low Resolution Imaging Spectrometer (LRIS) on 
the Keck I telescope and obtained a second spectrum on UT 2001 February 23.66. On 
both nights we used a 300 lines/mm grating blazed at 5000 A with a 1".5 wide slit. The 
effective spectral resolution FWHM varies from « 13 A on the blue side to ~ 11 A on the 
red side. A total of 1800 s of exposure were obtained on each night. The data were trimmed, 
bias-substracted and flat-fielded using standard procedures. Extraction of the spectra was 
performed using the IRAF APALL package. Telluric lines were corrected by fitting the 
continuum including the A and B atmospheric bands of a spectrophotometric standard. The 
spectra of spectrophotometric standards Feige 67, Hiltner 600, and HZ 44 were used for flux 
calibration and a Hg-Kr-Ne-Ar-Xe comparison lamp provided the wavelength calibration. 
The wavelength cahbrations have a typical accuracy ofRiO.2 — O.SAas determined by the 
root-mean-square (rms) deviation of the individual lines from the smooth function used to 
flt them. The large range of precision reflects the presence of fewer arc lines on the blue side. 
Finally, the LRIS-B (blue channel) and LRIS-R (red channel) spectra were connected at the 
crossover wavelength of the dichroic beam splitter. 

2.2. ESI Spectroscopy 

A higher-resolution spectrum of GRB 010222 was obtained on UT 2001 Feb 23.61 using 
the Echellette Spectrograph and Imager (ESI) mounted on the Keck II telescope. This mode 
spans ten orders with effective wavelength coverage from 3900 A to 10900 A. The spectral 
resolution is 11.4 km s^^ pix^^ with typical dispersion in the instrument ranging from 0.16 
A pixel"^ in order 15 to 0.30 A pixel^^ in order 6. In the echellette mode the instrument 
provides a wavelength resolution R 10000 over the entire spectrum. A total of 3 x 1200 
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second exposures were obtained and reduced using Tom Barlow's Mauna Kea Echelle Ex- 
traction program (MAKEE)^. MAKEE, originally designed for the High Resolution Echelle 
Spectrometer (HIRES), is an automated procedure that allows standard processing as well 
as wavelength and flux normalization. For the wavelength calibration wc used a CuAr cath- 
ode lamp which provides a large number of unsaturated lines to flt the spectrum with rms 
~ 0.09 A. The individual spectra at a common epoch were average- combined after running 
through the pipeline extraction. 



3. Reductions 

3.1. Absorption System Identifications 

We have identified three absorption systems (Bloom et al. 2001; Castro et al. 2001a) 
which were reported independently by Jha et al. (2001). The systems are located at 
zi — 1.47688, Z2 — 1.15628, and zs — 0.92747, where the redshift reported is the average of 
individual lines for each system. The line identification corresponds to a cross-correlation of 
features in the spectrum with a listing of known absorption lines (Verner, Barthel, & Tytler 
1994) assuming possible redshift systems at Zi, Z2, and z^. In order to weed out spurious fea- 
tures possibly due to sky subtraction, only those lines present in three spectra are reported. 
Moreover, we demanded Arest[l-0 -t- z(l-a;)] < Xobs < Kest[^-^ + z(l+a;)] and the correspond- 
ing oscillator strength ratios for identified fines. In the previous, x stands for a tolerance 
factor in the line identification and accounts for uncertainties in the wavelength calibration. 
Figure 1 shows the LRIS spectra including line identification for each absorption system. 
The higher-resolution ESI spectrum with identified absorption features is shown in Figure 
2. In the highest redshift system we find kinematic evidence for two distinct components 
at 2;ia=l. 47590 and zif,=l. 47688 respectively, with a rest-frame separation of ~ 119 km s^^. 
The kinematic structure of this system will be discussed at greater length in §4.3. Tables 1 
and 2 summarize the line identifications for each setup and list the measured vacuum wave- 
lengths, observed wavelengths converted into vaccum and corrected for heliocentric motion, 
estimated redshift, oscillator strengths fij, equivalent widths {Wo) in the rest frame for both 
nights, and error estimates in their equivalent widths. We have adopted the compilation of 
wavelengths and oscillator strengths by Morton (1991) with the revisions proposed by Savage 
& Scmbach (1996). The main source of uncertainty in the equivalent width measurements 
derives from the identification of the continuum level. In order to compute the errors in the 
equivalent width for each line we used the IRAF SPLOT task which allows error estimates 
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based on a Poisson model for the noise. For blended lines, IRAF SPLOT fits and deblends 
each line separately using predetermined line profiles. Error estimates for blended lines are 
computed directly in SPLOT by running a number of Monte-Carlo simulations based on 
preset instrumental parameters. 



3.2. Continuum Fitting 

The flux-calibrated spectra were corrected for Galactic foreground extinction assum- 
ing E{B — V) = 0.023 (Schlegel, Finkbeiner, & Davis 1998) and the extinction curve of 
Cardelli, Clayton & Mathis (1989). After correcting for Galactic extinction, the spectra 
were converted to F^, and u units and the continuum divided in equal size bins, excluding 
regions containing absorption lines. The statistical la error was determined by estimating 
the standard deviation within each bin. Next we found that a power law of the form Fj, oc z/'^ 
produced an adequate fit to the optical spectrum. The best fitted model changed slightly 
from (3 = -0.89 ± 0.03 on Feb. 22.66 to /3 = -1.02 ± 0.08 on Feb. 23.66. The contin- 
uum points along with the best fitted models are shown in Figure 3. Despite the apparent 
steepening of the spectra, there is one important caveat, the difficulty of doing spectropho- 
tometry through a narrow slit. Because we have used a 1.5" slit width the throughput of 
the spectrograph is reduced. Spectrophotometry might also be affected by atmospheric dif- 
ferential refraction (Filippenko 1982). On Feb 22.66 this effect was minimized by aligning 
the slit at the parallactic angle for the GRB as well as for the spectrophotometric standards. 
Moreover all the objects were observed at similar airmasses. During Feb 23.66 only one slit 
angle was used throughout the night and the airmasses for the spectrophotometric standards 
were slightly higher at the time of acquisition. Therefore, the uncertainty in the spectropho- 
tometry during the second night still allows the possibility that the spectral index did not 
change between the two nights. The difficulty of achieving accurate spectrophotometry in 
this case gives more weight to the results obtained from simultaneous multi-band photom- 
etry. Our derived spectra index f5 = —0.89 ± 0.03 appears in agreement with the reported 
single power law (3 = —0.88 ± 0.03 obtained from BVRI optical data (Stanek et al. 2001). 
Moreover the spectral behavior of the lightcurves does not show any steepening of the index 
P as a, function of time (Stanek et al. 2001). The previous result is also consistent with 
P — —0.89 ±0.03 derived from a low-resolution spectrum obtained five hours after the burst 
(Jha et al. 2001). But our spectral index is slightly shallower than the index (3 — —1.1 ± 0.1 
reported by Masetti et al. (2001). One possible cause for the steepening in the Masetti et al. 
(2001) results is the inclusion of the f/band data in the overall fit (Stanek et al. 2001). The 
bulk of the evidence seems to argue for a single power law spectral index for the duration of 
the burst. Hereafter we adopt (3 = —0.89 ± 0.03. 
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Our optical spectral index is similar within uncertainties to the measured X-ray spectral 
index — 0.97± 0.05 reported for the whole BeppoSAX observation (in't Zand et al. 2001). 
However, as pointed out by these authors, the optical flux is too faint or the X-ray flux is 
too bright to accommodate a simple synchrotron power-law. One possible explanation for 
this mismatch is an X-ray excess above the standard synchrotron spectrum due to inverse 
Compton (IC) scattering (Sari & Esin 2001). Harrison et al. (2001) have recently presented 
broad-band observations of GRB 000926 that are consistent with this IC interpretation. An 
alternative explanation for the difference between the optical and X-ray flux is reddening of 
the optical spectrum. The steep slope of the continuum fit to the optical spectra in addition 
to the absence of a 2175 A "bump" ubiquitous in a Galactic extinction curve implies that 
the extinction curve intrinsic to GRB 010222 would have to be approximately flat in order 
to not redden the spectra dramatically. This in turn suggests that the dust content at the 
burst location is relatively low or that any dust present mimics the gray variety (Aguirre 
1999). IR and submillimeter observations of GRB 010222 by Frail et al. (2001b) concluded 
that the majority of the reprocessed emission seen at the host galaxy is from a separate 
star burst component. By their estimation thermal radiation or dust scattering near the 
burst site cannot reproduce the bulk of the observed sub-millimeter energy. Gonsidering the 
possibility that the host is a starburst galaxy, we estimated the intrinsic reddening using 
the extinction curve of Galzetti et al. (2000). Dust corrections are particularly dangerous 
given the complexity of the dust geometry in any particular line of sight (Witt et al. 1992). 
Nonetheless, even a small color excess, E{B — V)^ 0.1, requires a flatter intrinsic power law 
for the afterglow with j3 = —0.39. This color excess is required to satisfy (3 = —{p — l)/2 for 
z/ < i/g in the optical band and /? = — p/2 at frequencies u > in the X-ray regime. Here p 
is the index of the power-law electron energy distribution, and Uc is the "cooling frequency" 
at which the electron energy loss time scale is equal to the age of the shock. The spectrum, 
corrected for such extinction at the host galaxy, is illustrated in Figure 4. Alternatively, it 
is possible that i/c has moved below the optical band at the time of our first spectrum, 8.45 
hr after the burst. If this is the case, the reddening argument does not apply; instead, an 
IG scattering interpretation is possible. 



4. Analysis 

4.1. Metallic Absorption Systems 

The statistical properties of metallic absorption systems have been studied extensively 
at various spectral resolutions {e.g., Steidel & Sargent 1992; Ghurchill et al. 1999). These 
studies have provided insight into the equivalent width distribution, evolution, and cosmo- 
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logical clustering of absorbers. In the particular case of GRB 010222, the equivalent widths 
of the Mg II doublet in the identified systems indicate that they are unusually strong in 
comparison to the population of absorbers studied by Steidel & Sargent (1992). The Mg II 
doublet ratios l^o(Mg II 2796.35)/l^o(Mg II 2803.53) are 1.06 ± 0.06, 1.17 ± 0.19, and 1.46 
± 0.12 for the 2;i=l. 47688, 2:2=1.15628 and 2:3=0.92747 systems, respectively. In each case 
the ratio differing from 2 implies that the Mg II lines are strongly saturated. In addition to 
the saturation of Mg II lines, we shall show evidence for saturation effects in the majority of 
Fe II lines in the ESI spectrum. We avoided using these hues for direct column density mea- 
surements given the likelihood of a non-standard curve of growth under saturated conditions 
{e.g., Morton & Bhavsar 1979). 

Instead, resolved absorption lines in the ESI spectrum that are not strongly saturated 
are more appropriate for column density determinations. The Mn II triplet at 2576.88 A, 
2594.50 A, and 2606.46 A in the 2;i=1.47688 system meets this requirement and appears 
particularly well suited to studying the highest column absorption system at zu in which it 
resides. These lines fall in the region of intermediate optical depth, r > 1, i.e., on the "flat" 
part of the curve of growth. Here the observed broadening of the Mn II lines is assumed to 
be due to the random motion of the gas, described by a Maxwellian velocity distribution. In 
this regime the equivalent width becomes sensitive to the Doppler parameter b, which can 
be determined directly from the broadening of the absorption lines. The measured Doppler 
parameters b are listed in Table 3. In order to derive column density Nj for any metal j, we 
used the general prescription for a standard curve of growth (Spitzer 1978) 



where Nj is written as a function of equivalent width Wx, rest wavelength A, oscillator 
strength fij, optical depth r, and its dependence on the Maxwellian velocity distribution 
F(r). Metallic column densities log(A''j) as well as equivalent neutral hydrogen column 
densities log(A''Hi) assuming solar abundances (Anders & Grevessc 1989) arc given in Table 
3. Notice the overall agreement in the column densities derived using three different Mn II 
lines. The column densities listed in the fifth and sixth columns of Table 3 have not been 
corrected for depletion into dust grains (Jenkins 1987), which is an uncertain quantity. 

Waxman & Draine (2000) have shown that dust may be sublimated by a prompt opti- 
cal/UV flash at any distance closer than fa 10 pc from the burst. Fruchter, Krolik & Rhoads 
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(2001) also considered the possibility of grain charging as a mechanism for dust destruction. 
This could take place as far away as ~ 100 pc from the birthsite. Thus, the combined effect 
of both destruction mechanisms, temperature of the medium, and the return of dust depleted 
metals to the gas phase make the depletion correction largely uncertain. For completeness 
we make the simplest assumption ignoring any immediate return of metals into the gas phase 
and adopt a depletion correction of —1.45 dex for Mn"'' based on Goddard High- Resolution 
Spectrograph (GHRS) measurements towards ( Oph obtained by Savage & Sembach (1996). 
The resulting corrected values are shown in the last column of Table 3. Given the uncertainty 
in the amount of depletion it appears safe to regard the uncorrected HI column densities 
as lower hmits. We note that the weighted average of Mn+ corrected for depletion implies 
log(A^Hi)corr = 21.89 ±0.13 cm~^ in the Zia system, which is similar to the rcdshift-corrected 
column density log(A''Hi) ~ 22.00 cm~^ derived by in 't Zand ct al. (2001) from the X-ray 
afterglow. The results presented so far assume solar abundances, alternatively, less than 
solar abundances for Mn 11 might accommodate the difference between metal abundances 
and the values reported by in 't Zand ct al. (2001) and Salamanca et al. (2001). 

We have completed the estimates of column densities Nj by including blended absorption 
lines. The term "blended" is used here to refer to the lines where the Zia and Zn, systems 
cannot be separated cleanly in the ESI spectrum. In particular, by comparing the equivalent 
widths and the oscillator strength ratios for Zn II (2026.14 A, 2062.66 A) and Cr II (2056.25 
A, 2062.23 A, 2066.16 A) lines we see that they fall on the linear portion of the curve of 
growth. For these lines we can obtain a direct measurement of the column density. On 
the other hand, the Fe+ and Mg+ hues, although structured, are strongly saturated. For 
that reason we decided to fit Gaussian profiles at the positions of the identified zia and zu 
components while running the deblending routines. In this case the linear part of the curve 
of growth only provides a lower limit to the column density Nj through the equation 

.3 xlO'^gg ,3) 

The resulting column densities are listed in Tables 3 and 4. In each case the density corre- 
sponds to the average of various single transitions. As a consistency check, the HI column 
density derived assuming solar abundances for the generally undepleted Zn+ is in agreement 
with the values cited in Table 3 for the resolved Mn+ lines. 
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4.2. Is ^1=1.47688 the redshift of GRB 010222? 

The absence of any emission lines in the spectrum shghtly compUcates the unequivocal 
identification of this system with the host galaxy. However, there is information in the 
optical spectrum that supports this interpretation. First, the non-detection of a Lyman-a 
break places an upper Hmit for the GRB host galaxy of z < 1.96. Second, Salamanca et 
al. (2001) have reported evidence of a red wing of Ly-a absorption at z 1.476. Third, 
the measured equivalent widths of the Mg II doublet and the ratio of Mg I to Mg II in this 
system are also among the highest known for GRB afterglows and for the sample of Mg 
II absorbers compiled by Steidel & Sargent (1992) at similar spectral resolution. Fourth, 
we can quantify the probability of detecting three intervening galaxies in the line of sight. 
The LRIS setup used during the first two nights spans the 3600-9205 A wavelength range 
with a rest equivalent width threshold limit Wo > 0.3 A, which matches the resolution and 
equivalent width detection limit of the Steidel & Sargent survey. This wavelength coverage 
can be translated into a redshift range for Mg II lines of 0.284 < z < 2.284, or a redshift 
interval Az — 2.0. Thus, the number of absorbers per unit redshift is N/Az — 1.5, which is 
less stringent than an earher estimate of this quantity (Jha et al. 2001), but is still larger 
than the expected number of Mg II absorbers per unit redshift (N/z) = 0.97 ± 0.10 for a rest 
equivalent width threshold of 0.3 A (Steidel & Sargent 1992). The probability for the chance 
superposition of three absorbers in this line of sight is non-negligible but small. Lastly, if the 
highest redshift of an Mgll system is 2;= 1.47688, then, in 80% of the cases, the source could 
have originated from no further than z ~ 1.8. If it originated from a larger distance then, 
given the rapid increase in the number of metal absorption systems with increasing redshift, 
a source redshift higher than z=1.8 would be unlikely (Bloom et al. 1997). Although not 
definitive proof, these arguments are most consistent with the 2;i=l. 47688 absorption system 
as the host galaxy of GRB 010222. 



4.3. Kinematics and abundances of the 2;i=l. 47688 absorption system 

The advantage of using higher spectral resolution is the capability to resolve multiple 
components along the line of sight. In this instance, we identify two kinematic components 
at 2;ia=1.47590 and ^16=1. 47688 in the ESI spectrum. The presence of multiple components 
immediately raises the question of structure within the host galaxy. We might be probing a 
galactic disk accompanied by an additional kinematic component introduced by absorbing 
clouds in the halo (Disk-Halo scenario), or instead the two components could be due to 
absorption arising in distinct regions of a galactic disk (Disk scenario). In order to test 
these scenarios we also used information about element abundances and dust depletion. In 
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particular, the ratios [Zn/Cr] and [Zn/Mn] provide a reasonable measure of the gas-to-dust 
ratio and the temperature of the absorbing medium. 

The main hne of argument lending support to a Disk-Halo scenario is the absence of 
a Zia kinematic component in the resolved set of Mn+ lines. This is better illustrated in 
Figure 5, which shows the kinematic components for different metallic lines in velocity space 
using zib = 1.47688, which is assumed to be the systemic redshift, as the velocity zero point. 
Mg°, Mn+ and Mg+, are present in the Zu, system, while the Zia component lacks Mn"*" but 
contains Mg^ and some of the structured but likely saturated Fc^ lines. Since the ionization 
potential of Mn"*" is almost identical to that of Mg"*", a possible reason for the absence of Mn"*" 
at Zia is an underabundance of Mn in a metal-poor environment like the halo. This is based 
on the conjecture that Type la supernovae are an important source of Mn (Samland 1998, 
Nakamura et al. 1999), or that there is a metallicity dependence yield of Mn in massive 
stars (Timmes, Woosley, & Weaver 1995), as discussed in the case of damped Ly-a systems 
by Pettini et al. (2000). An early appearance of type la SNe during the halo phase in Milky 
Way formation models has been questioned (Chiappini, Mateucci, & Romano 2001). If the 
latter authors are correct, our Zia component might arise in a halo cloud, separated from the 
velocity zero point by ~ —119 km s~^, with signatures of Mg+ and Fe+ but underabundant 
in Mn+. 

The possibility of a Disk-Halo line of sight makes the abundance ratios for each kinematic 
component that much more crucial in probing the ISM at the host galaxy. Zinc is a fairly 
reliable metallicity indicator since it resides more than any other metal in the gas phase. 
On the other hand, metals like Cr and Mn are more easily depleted onto grains and a large 
fraction of their abundance is contained in metallic dust grains in the galaxy. If the Zib 
component is an absorption feature of the halo then a higher gas-to-dust ratio than in the 
disk is expected. Unfortunately, the [Zn/Mn] ratio in the zu component is only an upper 
limit given that Zn+ lines from both velocity components are blended, as are the Cr+ lines. 
We note however that the Mn and Cr abundances are comparable. In fact the observed 
[Zn/Cr] ~ 0.61 is less than 1.61 > [Zn/Cr] > 1.04 observed towards ( Oph (Savage & 
Sembach 1996) and 1.85 > [Zn/Cr] > 0.75 measured in a sample of 20 sightlines observed 
with HST (Roth & Blades 1995). Although the [Zn/Cr] ratio is less than in any measured 
Galactic sighthne it does not represent a clear discriminant between a halo and a disk or 
a damped Lyman absorber since the spread in [Zn/Cr] is large for either sample (Roth & 
Blades 1995). 

Our results might be an indication of low dust content but not at the extreme reported 
by Salamanca et al. (2001). It seems that their Zn II measurements are systematically 
lower than those presented here or any other pubhshed results for GRB 010222 (see Jha 
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ct al. 2001, Masetti et al. 2001). In addition, Zn 11(2026.14 A) listed at Xobs = 5011.2 
A in Table 2 of Salamanca ct al. (2001) appears to be a misidcntification. We find that 
line is more consistent with Fe 11(2600.17 A) at the system while Zn II should be around 
^obs = 5018.98 A at the higher redshift. For completeness we mention that we have not 
detected any of the unidentified lines listed in Table 3 of Salamanca et al. (2001) in the 
region covered by our spectra. Finally, the reduced [Zn/Cr] ratio in the spectrum might be 
hinting that we have detected a warm medium in the zij, component where depletion for Cr 
and Mn converge (Savage & Sembach 1996) or else that there is httle dust to deplete onto 
in this environment. A warmer medium would result in different dust depletion corrections 
and lower the inferred values of A^hi assuming standard solar abundances. These results are 
in agreement with the depletion patterns discussed by Savagho et al. 2002. 

We now consider the alternative of a Disk scenario. In general, any neutral Mg present 
in a galactic disk must be shielded from external photoionizing flux. A feasible scenario 
is that we are probing a line of sight passing through two distinct regions within the disk 
containing neutral gas. These regions could be underdepleted molecular clouds cocooning 
Mg in its neutral form. A possible consequence of uneven photoionization in the disk is 
the marginal detection of a blue "wing" on the Zia component of the Mg II (2796.35 A) 
absorption line (Figure 5). This feature would argue for a stronger photoionizing source 
near the closer side of the disk and against a galactic halo. This option would eliminate the 
necessity for a halo component, but it does not explain the absence of Mn+ in the zia system, 
assuming similar enrichment processes for different parts of the disk. Instead, there may be 
a strong metallicity gradient along the disk, as has been observed along the disk of the Milky 
Way (van Steenbcrg & ShuU 1988). Assuming that a Disk scenario is correct and that the 
[Zn/Mn] ratio peaks towards the inner part of the host galaxy; the velocity systems might be 
tracing a fraction of galactic rotation. Independent of the interpretation, the presence of a 
galactic disk appears necessary to explain the properties of the strongest absorption system. 
This furthers the connection of long-duration GRBs with a disk progenitor. 



4.4. Looking for signatures of H2 absorption and Fluorescence 

The impact of the associated UV fiash on H2 molecules in the vicinity of a GRB has 
been determined by Drainc (2000). The direct observable consequence of strong UV radiation 
from a GRB on the surrounding H2 is the production of vibrationally excited levels which 
could create strong line absorption in the region 912 A < Arest ^ 1650 A, and reradiatcd 
fluorescent emission in a similar range. For a redshift of z = 1.477 this corresponds to the 
region Aobs < 4087 A. Unfortunately the region where the absorption effect is strongest at 



- 12 - 



Aobs ~ 3096 A falls outside our spectral coverage. Thus, we are only able to probe a fraction 
of the absorption region near the blue end of the spectrum where the quantum efficiency 
is low and the noise is high. The strength of the absorption depends on the amount of H2 
available for vibrational excitation. At the LRIS resolution we find no definite evidence of H2 
absorption for Aobs < 4087 A to make a proper determination of Nh2- However, our spectra 
probably rule out excited states of H2 at > 10^° cm^^ based on the transmission plots 
of Draine (2000). 



4.5. Absorption-line variability 

If the line absorption occurs sufficiently close to the burst site, the afterglow emission 
could manifest itself through photoionization of the GRB environment, thus, variability of the 
absorption features. To detect this time-dependent effect various snapshots of the absorption 
spectrum are required. In this particular instance we take advantage of the low-resolution 
spectra obtained over two consecutive nights to search for the predicted effect. The measured 
equivalent widths for each night as well as the statistical uncertainties are listed in Table 
1. Comparison between the LRIS equivalent widths on consecutive nights shows no strong 
evidence for time dependence of the absorption lines within statistical uncertainties. There 
is also an overall agreement with the ESI measurements. We deduce that discrepancies in 
any two measurements arc mainly due to uncertainties in the continuum level in the vicinity 
of each line. There does seem to be a systematic trend of larger equivalent widths on the 
second night below 5403.80 A. However, we also note that it is clear from the spectrum that 
the equivalent widths of lines on the blue side are relatively "weak" with respect to well 
detected lines on the red side where such trend is not detected. Hence, at this point we 
cannot attribute the systematic trend entirely to a physical effect since the relative line flux 
measurements might have been affected by the continuum fitting during the second night as 
discussed in §3.1 and §3.2. 

Although the equivalent width depends on the overall characteristics of the continuum, 
its main uncertainty is defined locally by the error in determining the proper level of the 
continuum on either side of the line of interest. As discussed in §3.1 the errors quoted for 
equivalent widths correspond to a Poisson model of the noise in the region of interest. Such 
an assumption might be oversimplified and not account for true variations in the continuum, 
but wc believe that wc safeguard against strong systematic effects by using the same, con- 
sistent method in all measurements. Hence based on our analysis we do not find statistically 
significant changes in our measurements. Taking this comparison further there are pub- 
lished equivalent width measurements for this burst elsewhere (Jha et al. 2001, Masetti et al 
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2001, Salamanca et al. 2001). Although the duration of observation, spectral resolution and 
brightness of the burst varies in each case, the reported measurement techniques are similar. 
Direct comparison finds that our results are consistent with previous measurements. Never- 
theless some differences arc present in a few cases such as C IV, Cr II and Zn II that were 
discussed in §4.3. However, wc do not find systematic and statistically significant changes 
as a function of time in the published data that should be present if either photoionization 
or dust destruction in the circumburst medium are actively dominant for similar species. 
We refrain from reporting a measurement for Al 111(1862.79 A) during the second LRIS 
spectrum in Table 1 because of the reduced signal-to-noise in that region of the spectrum. 
For completeness we note that the hne has not disappeared and it is clearly detected in the 
ESI spectrum obtained at a comparable time (Table 2) . 



In order to interpret our observations under contraints imposed by the lack of absorption- 
line variability, we simulated the GRB environment using a standard photoionization model. 
The photoionization evolution has been described by Perna & Loeb (1998) and Bottcher et 
al. (1999). The photoionization rate can be written in the following form: 



where rij is the density for species j, F^/{r,t') is the flux at radius r, t' is the elapsed time 
in the GRB frame and crbfj(i^') is the photoionization cross section at frequency u' using the 
analytical fits provided by Verner et al. (1996). We neglect recombination processes since 
the densities we will be considering are typically n < 10^ cm~^, where the corresponding 
recombination timcscalcs are much longer than the duration of the GRB and our observations 
(Perna & Loeb 1998). In addition the full effect of dust destruction and the return of depleted 
metals into the gas phase has been ignored in the code. 

Next we constructed the input flux F,y! (r, t') for this particular afterglow from the ob- 
served time decay (Stanck et al. 2001) and the measured optical spectral index. The function 
has two forms to accommodate the observed break in the light-curve at tbreak=0.72 day. Thus 

for t < 0.72 day 



5. 



Photoionization Code 
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F,,{ro,t') = 5.26x10-28 ( — ^ ( -^-^) ( i^Il±^^ ''"ergs cm-^ sec"! Hz"! 

^ ^ V6-94 X 1014(1 + ^)Hzy V(l + '2)^0/ V 0-352 day/ ^ 

(5) 

otherwise 

F,iro,t') = 2.97x10-28 f — ^ f ■^^^) ( ^^il±^^ ''"ergs cm-^ sec-^ Rz^' 

^ ^ V6.94 X 1014(1 + z)Hzy V(l + -2)^0/ V 0-72 day/ ^ 

(6) 

where c? is the luminosity distance to the burst at z = 1.47688 (assuming Ho = 65 km s^^ Mpe-i, 
flm = 0.3, A = 0.7), and Tq = 10^^ cm is the inner radius of the photoionized region set 
as a boundary condition. We have used the observed flux and ignored any intrinsic dust 
corrections. In a previous section we showed that intrinsic reddening following a Calzetti 
et al. (2000) extinction law for starbursts would dim and stepeen an intrinsically fiat spec- 
trum. The effect of using an input flux uncorrected for intrinsic extinction is to slow down 
the time-dependence in our models; therefore, the integration times presented here should 
be considered upper limits. The input flux has been integrated in the UV range up to a 
cutoff energy of 0.2 keV. The reasoning behind this choice is partly arbitrary given the fail- 
ure to locate a precise spectral break at the "cooling frequency" Pc between the optical and 
X-ray ranges. The selected cutoff energy of 0.2 keV represents roughly the region where the 
photoionization cross section for the species of interest has dropped sharply. Above 0.5 keV, 
X-rays are mostly absorbed by heavy elements. 

Wc start the photoionization of a pristine neutral medium at initial time t'^ that follows 
the detection of prompt emission in the BcppoSAX energy band. The transition between 
the prompt GRB emission and the onset of the afterglow is not a very well defined point 
observationally. Variability arguments set the scale for internal shocks to take place in a 
region of size R^i 3x lO^^ cm {6t/l s)(r/100)2. External shocks responsible for the afterglow 
become significant at ~ 10^^ cm from the source (Piran 1999). For numerical reasons we 
have selected an initial time of t'^ — 10 seconds for the onset of photoionization. 

The photoionization process can be pictured as ocurring inside a sphere of outer radius 
fmax: Subdivided into thin spherical shells of width Ar so that each shell i is optically thin. 
In each shell the initial incident flux spectrum described above will be absorbed and will also 
follow a r"^ law as r — > rmax- Hence the incident spectrum at any shell i will be given by: 
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We ignore bound-bound processes therefore r,^,, stands for the photoionization optical depth 
which is estimated within each shell i as: 



where r,^^j << 1. For our models we included H, He, and Mg and neglected other elements. 
The initial number density rij for each species is largely unknown but we decided to use 
standard solar-system abundances. For simplicity we also assumed that initially Mg is equally 
divided into Mg° and Mg+^ 



Since the depletion correction/abundance and therefore the total A^h near the burst 
site is not well determined, we examined the time dependence of equivalent widths for two 
different hydrogen column densities. The first, log(A^Hi) = 20.44 cm~^ represents the inferred 
HI column ignoring any dust depletion correction and takes into account the possibility of 
complete circumburst dust destruction. The second, log(iVHi) = 22.00 cm^^ is appropriate 
for standard Milky Way depletion onto dust or perhaps less than solar metal abundances. 
The latter is also consistent with measurements of the red wing of a possible Lja feature 
(Salamanca et al. 2001), and the BeppoSAX afterglow results (in't Zand et al. 2001). In 
each case we chose four different combinations of number densities rim and rmax to reproduce 
the simulated hydrogen column density A^hi = nm x rmax- The parameters for each model 
are listed in Table 5. The denser, more compact models represent the conditions in dense 
molecular clouds which have been postulated as a potential ambient environment for the 
progenitors of GRBs. On the other hand, the low density and extended environments are 
compatible with low density molecular clouds, young stellar clusters, and possibly galactic 
superbubbles (Scalo & Wheeler 2001). 

We considered it relevant to trace the time evolution of two particular absorption lines 
present in our spectra, Mg I (2852.96 A) and Mg II (2796.35 A). Mg I is an appealing line 
given the possibility of encountering Mg in its neutral form in regions of the galaxy shielded 
from photoionization. Mg 11(2796.35 A) is of interest because of its detection in a number of 
GRB spectra. An additional observational reason for choosing these particular lines to test 




(8) 



a 



6. Numerical Results vs. Physical Conditions 
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our models is the consistency of the hne flux measurements on the red side of the spectrum 
where the effects of differential refraction appear of less importance. In order to model our 
observations, we integrated to a maximum time tmax = 32.45 hours after the GRB, when the 
second LRIS spectrum was obtained. The equivalent width at different stages was estimated 
using equations (1) and (2) which return Wx. Throughout our calculations we have assumed 
a constant Doppler parameter b representative of the system (see Table 3). This appears to 
be a fairly robust assumption since the time dependence for the lines is mainly a function of 
the photoionization equations. 

Figure 6 illustrates the time evolution of the equivalent width for Mg 1 (2852.96 A) 
and Mg 11 (2796.35 A) in a log(iVHi) = 20.44 cm~^ environment. The times when the LRIS 
spectra were obtained are indicated. Notice the distinct differences in the line evolution in 
each case. The denser, more compact models are ionized rather rapidly since most of the 
gas is close to the GRB and gets impacted by a large flux. On the other hand, less dense 
models representative of thin galactic environments require longer for noticeable effects to 
take place. The nm — 89 cm~^, rmax = 1 pc model appears to mark the line where our low- 
resolution spectroscopy would have been able to detect time dependence with this column 
density. Denser environments than this appear to be disfavored according to our model. 

Figure 7 shows the time dependence for a log(iVHi) = 22.00 cm~^ environment. Similarly 
compact, high density environments are photoionized rather quickly while extended, less 
dense ones would evolve with a shallower decline. The n^i = 2 x 10^ cm^^, rmax = 5.0 x 10^^ 
cm model is close to our detection limit for this denser environment. The results in this 
case are more dependent on the choice of the inner radius Vg because the scale of some 
of the models becomes comparable to Tq. A different choice for the inner radius would 
either accelerate or delay the photoionization timescale, therefore, compact models with 
this column density are more tentative. However, it appears that very dense, compact 
environments would have been more highly ionized than is observed in our initial spectrum. 

One interpretation of X-ray and optical afterglow data from GRB 010222 invokes a 
transition from a relativistic to a non-relativistic regime, and consequently the deceleration 
of the shock, to account for the break in the light curve (in 't Zand et al. 2001). This 
scenario requires an environment with nm ~ 10^ cm~^ and size ~ 10^^ cm. Although 
our simulations argue against the presence of a dense 10^^ cm neutral gas slab between 
the two observations, we cannot rule out the presence of a large density slab at the time 
of the burst itself. However our modelling also presents an alternative scenario for more 
extended absorption in the host galaxy rather than a single slab. Recent modelling of 
the multiwavelength observations of GRB 010222, including the radio emission, finds that 
a jet model expanding into homogeneous external medium with density nm ~ 3.2 cm~^ 
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is consistent with the broadband data (Panaitescu & Kumar 2001). A jet and low-density 
interpretation is also required to satisfy the standard-energy output with beaming corrections 
of Frail et al. (2001a). One final point in this discussion is the shape of the 7-ray burst light 
curve and baryon contamination. A very dense nearby medium gives rise to high optical 
depth r where baryons could prevent relativistic expansion. The overall effect could lead to 
"dirty" fireballs and smearing of the 7-ray profile. Instead the burst profile appears highly 
structured (in 't Zand et al. 2001). 

Consider now the actual column density to this burst if the host galaxy has a dust- 
to-gas ratio similar to that of the Milky Way. A color excess of E{B — V) = 2.00 mag 
is estimated for log(iVHi) = 22.00 cm^^ and assuming a Galactic conversion Niii/E{B — 
V) — 5.0 X 10^^ cm"^ mag~^ (Savage & Mathis 1979). This is significantly larger than 
E{B -V) ^ 0.06 mag for a log(A^Hi) = 20.44 cm"^ environment. It is possible that the UV 
emission from the burst is responsible for destruction of neighboring dust in the high-density 
case (Galama & Wijers 2001), but further extinction at larger distances is expected along 
the line of sight unless the burst is near the Earth-facing side of the host galaxy (Reichart & 
Price 2002). We must keep in mind, however, the possibility that we are probing different 
scales in two separate energy bands. 

7. Implications for the Nature of the GRB 010222 Progenitor Environment 

GRB 010222 is a good candidate to explore whether we can place constrains on the 
progenitor environment. The kinematic components resolved in the high-resolution spectrum 
give support to a location of this burst in a galactic disk. High-resolution spectroscopy of 
GRB 000926 (Castro et al. 2001b) favors a similar interpretation for that burst. The slope 
of the spectral continuum in this instance suggests that the burst occurred in a dust-free 
environment or was surrounded by a screen of gray dust (Aguirre 1999). Near-IR and 
submiUimeter observations by Frail et al. (2001b) also appear to rule out the possibility 
of strong dust scattering or thermal emission near the GRB site. However it is difficult to 
model a priori the properties of dust at this redshift since the distribution and properties of 
galactic dust as a function of redshift are largely unknown. 

7.1. Isolated Molecular Clouds 

Popular models of gamma-ray burst progenitors include young massive stars (Woosley 
1993, Vietri & Stella, 1998). This makes the interior of molecular clouds a potential birth- 
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place for GRBs. Molecular clouds tend to comprise large amounts of dust and clumpy 
pockets of gas ideal for star formation. If GRB 010222 took place in a molecular cloud, the 
slope of the optical continuum indicates that either the nearby dust was destroyed rather 
immediately by the burst through a combination of sublimation by an optical/UV flash 
(Waxman & Draine 2000) and grain charging (Fruchter, Krolik & Rhoads 2001) or that 
the surrounding dust content was low prior to the explosion. The absence of measurable 
signatures of metals being returned to their gas phase as a by-product of dust destruction 
on such a short timescale does not allow direct distinction between these two possibilities. 
However, a high fraction of dust results in a larger H2 formation rate since molecular hy- 
drogen tends to form in reactions taking place on grains. The absence of strong absorption 
features due to vibrational excited levels of H2 expected from strong UV emission impacting 
the surrounding H2 that might permeate a progenitor in a dusty molecular cloud (Draine 
2000) is conspicuous. Moreover, assuming a mean density of 10^ cm~^ and the typical size 
of a dense molecular cloud, the absorption lines detected should have weakened over the 
period of our observations. A massive and dense molecular cloud might also give rise to 
reflection at 6.4 keV as seen near the Galactic center (Murakami, Koyama, & Maeda 2001). 
Nevertheless, it is still viable that a molecular hydrogen detection has been made blueward 
of our spectral limit by Lee et al. (2001) and Galama et al. (2002). As reported by these 
authors the Atix deficit in the U band which extends below 3500 A does not fit the SMC or 
LMC-like extinction models. A depression signature is also expected if there is absorption by 
a Lyman-a wing (Salamanca et al. 2001). The absence of molecular hydrogen and reduced 
depletion resulting from the [Zn/Mn] and [Zn/Cr] ratios still allows the location of the GRB 
to be near the edge of a molecular cloud prior to the explosion. Thus, at the time of the 
event the fraction of gas to dust can be larger than at the center of the cloud, inhibiting 
the formation of H2. Immediately after the onset of the UV emission from the burst, any 
remaining dust could have been destroyed efficiently by the expanding UV radiation pulse, 
minimizing reddening of the optical continuum. 



7.2. Superbubbles 

Panaitescu & Kumar (2001) in their multiwavelength analysis of well-studied afterglows 
found that the ambient density for some bursts can be as low as 1.9 x 10~^ cm~^; however, 
some arguments favor higher ambient densities once the complete radio data set is included 
in the analysis of individual bursts (Berger et al. 2002). The presence of low densities led 
Scalo & Wheeler (2001) to argue for active star formation regions in superbubbles as possible 
environments for GRBs. Our photoionization models show that low-density environments 
with the typical sizes of superbubbles would not produce strong time- dependence in absorp- 
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tion lines. Moreover, the interior of a superbubble could be partionally ionized prior to the 
GRB event. The prior ionization might be accomplished via overlapping ionization bubbles 
generated by supernovae in clusters of massive stars. A pre-ionized environment decreases 
the amount of neutral material seen by the expanding photoionization front and prevents the 
full effect of time-dependent absorption. But it also enables the detection of time dependence 
in high-ionization species. We detect high-ionization C IV absorption in our low-resolution 
spectrum, but it is not clear that equivalent width variations are seen. In general C IV is 
a tracer of additional diffuse warm and hot gas that could originate in a superbubble envi- 
ronment. Unfortunately, our high-resolution spectrum does not cover this absorption line to 
firmly estabhsh if it comes from the disk of the host galaxy or rather from a gaseous halo. 
We note however that the equivalent width of C IV (1548.20 A) is comparable to that of the 
Mg II absorption lines. An additional signature of prior ionization could be the attenuation 
of the X-ray continuum by an ionized absorber like those seen in AGN spectra. The absence 
of strong absorption features except for a marginal Fe XIV-XVIII edge suggested by in 't 
Zand et al. (2001) prevents a straightforward connection from being made between X-ray 
and UV absorption features. 

An advantage of the superbubble scenario is the availability of a dust destruction mecha- 
nism and a favorable geometry for observation. The combination of winds and SN explosions 
would tend to blow the gas and dust away leaving filamentary or shell-like patterns like the 
one reported in the Galactic Orion- Eridanus bubble (Guo et al. 1995). In addition, winds 
and radiation pressure in the interior of a superbubble might selectively destroy small grains, 
a mechanism for gray dust production locally. However, this selective process is balanced by 
rapid destruction of large grains by the burst and earlier destruction by X-rays emitted inside 
the superbubble that might lead to grain charging (Fruchter, Krolik, & Rhoads 2001). In 
combination, these two processes could counteract the expulsion of small grains and suppress 
the formation of gray dust. Either dust destruction or its gray nature arc consistent with 
the slope of the observed continuum and the lack of strong dust reprocessing near the burst 
site (Frail et al. 2001b). The only practical problem with this interpretation is the mismatch 
between the optical and X-ray spectra. It is possible that the cooling frequency has moved 
out of the optical band. In this case we propose that the observed X-ray excess (in 't Zand 
et al. 2001) is mostly due to inverse Compton scattering rather than intrinsic reddening 
near the burst location but if true the ambient density cannot be below a few cm~'^ or in- 
verse Compton scattering would be suppressed (Harrison et al. 2001). If the effects of dust 
depletion are less pronounced than in the Milky Way, a log(A''Hi) = 20.44 cm~^ environment 
with rim < 89 cm~^ and rmax > 1 pc is suggested for GRB 010222. This is consistent with 
the best-fitted parameter nni ~ 3.2 cm~^ resulting from jet modelling (Panaitescu & Kumar 
2001). 
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7.3. Young Stellar Clusters 

The recently discovered Arches cluster (Nagata et al. 1995; Cotera et al. 1996; Serabyn, 
Shupe, & Figer 1998; Figer et al. 1999) consists of more that 100 young stars with masses 
greater than 20 Mq within a ~ 0.3 — 0.4 pc radius. Radio observations of the Arches 
stellar cluster have confirmed the presence of powerful ionized winds with mass-loss rates 
^ {1 — 20) X 10~^ Mq yr~^ coincident with massive stars in the cluster (Lang, Goss, & 
Rodriguez 2001). If a low metaUicity massive star that has lost its envelope turns out to 
be the prototype for the coUapsar model of GRBs (MacFadyen & Woosley 1999), a young 
stellar cluster environment appears to be a perfect cradle for progenitors as prelude to GRB 
formation. A compact young cluster environment is consistent with our observations and 
models but requires low dust content or rapid destruction near the burst site, as well as 
photoionization of the gas in the region by hot young stars prior to the GRB. The complexity 
of wind interaction and shocks might provide favorable dust viewing geometries and dust 
distribution similar to those discussed in the superbubble scenario. In addition, the UV 
absorption lines and their lack of time dependence might be attributed to the medium 
well outside the stellar cluster where the photoionizing flux is not strong enough to imprint 
substantial changes. One problem with this interpretation is the type of metallicities required 
by the MacFadyen & Woosley (1999) model which can be thought of as LMC-like (A.I. 
MacFadyen 2002, private communication). The abundances in the GRB 010222 host are 
larger than the abundances reported in LMC line of sights (Welty et al. 1999). In addition, 
the absence of strong signatures of ionized absorption and pronounced emission lines in the 
X-ray spectrum (in't Zand et al. 2001) might argue against a pre-ionized medium. 



7.4. Halo 

An alternative GRB progenitor is a coalescing compact object (Neutron Star-Neutron 
Star, Black Hole-Neutron Star) (Paczyhski 1986). Even though the distribution of birth kick 
velocities for compact objects has not been fully surveyed, typical velocities around a few 
hundred km s~^ could place a fraction of potential progenitors at high-latitude above the 
disk. The apparent presence of a galactic disk in this case seems to rule out the possibility 
of a halo origin for GRB 010222. Even though a location in the halo implies a low-density 
environment, the ionized regions of interest would be relatively compact clouds in the halo. 
In this case, it is difficult to reproduce the measured neutral column densities and observed 
time evolution in a halo event. Photoionization in halo events might be more effective in 
erasing absorption features entirely if the line-of-sight does not intercept the disk of the 
host galaxy. Current ideas suggest that the hard-short class of GRBs are due to coalescing 
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compact objects (Janka et al. 1999, Fryer et al. 1999). Rapid localization of the hard-short 
GRBs might allow the detection of absorption-line time dependence in high-ionization species 
following coalescence in a galactic halo. 

8. Conclusions and other Observational Reflections 

The use of high-resolution spectroscopy has resolved the host galaxy of GRB 010222 
into multiple components. We conclude that the GRB took place in a galactic disk that 

gives rise to the strongest absorption system. The power-law index of the optical continuum 
indicates that the dust content at the burst place is low or that the reddening follows a gray 
dust model. Under these conditions the bulk of the X-ray excess (in 't Zand et al. 2001) 
can be attributed to inverse Compton scattering. Low resolution spectroscopy obtained 
over the span of two days failed to reveal any significant evidence for time dependence of 
the observed absorption lines. Signatures of strong H2 absorption and/or fluorescence are 
also absent from our spectra. Photoionization models constructed to test a range of initial 
conditions for this burst show that dense, compact media would be photoionized fairly quickly 
while less dense media with large radii evolve slower with time. We argue that prior to a 
GRB, an environment partially ionized either by the progenitor or by nearby sources would 
prevent the full manifestation of predicted absorption-line evolution. The physical attributes 
of GRB 010222 are consistent with a low-density molecular cloud, superbubble, or young 
stellar cluster as the possible environment of the GRB. 

In the future some consideration should be given to a prompt optical-UV flash like the 
one accompanying GRB 990123. The intensity of such an event could in principle lead to 
dust destruction, but it might also be responsible for the photoionization of a fraction of the 
neutral gas. The full effect of such a flash has been ignored in our models. Preferentially, 
the observation of time-dependent absorption lines would work best in bursts occurring in 
the Earth-facing side of the host galaxy and in compact regions. Our analysis suggests that 
this technique applied with high-resolution spectroscopy could potentially distinguish the 
metallicity of progenitor environments. In that regard shortening the interval between burst 
localization and initial spectroscopy becomes crucial to exploring earlier, less photoionized 
times. The next generation of high-resolution imagers might directly resolve the environment 
of GRBs in a relatively nearby host galaxy. A prospect raised by the observations of GRB 
010222 is the study of the properties of galactic dust at high-rcdshift using bursts as beacons 
to illuminate the dust nearby. Further exploration of this issue is encouraged as spectroscopy 
is obtained for a larger sample of GRBs. 
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Table 1. LRIS Line Identification 



Line(Avac(A)) 


Ahelio(A) 


z 


fij 








Si 11(1526.71 A) 


3782.12 


1A773 


1.10x10' 


-1 


1.51 ± 0.13 


i. / Z 


1 

± 


U.4/ 


C IV(1548.20 A) 


3835.36 


1.4773 


1.91x10" 


1 

-± 


2.55 ± 0.23 


Z.D/ 


± 


u.zy 


+ C IV(1550.78 A)'' 






9.52x10" 


-A 










Fe 11(1608.45 A) 


3984.61 


1.4773 


6.19x10" 


9 
-z 


0.40 ± 0.04 


U.DZ 


1 


U.oU 


Al 11(1670.79 A) 


4138.77 


1.4771 


1.83 




1.48 ± 0.09 


1.58 


1 


O.lz 


Si 11(1808.01 A) 


4478.51 


1.4770 


5.53x10" 




0.60 ± 0.08 


U.oi 


_(_ 
It 


u. io 


Al 111(1854.72 A) 


4592.77 


1.4763 


5.60x10" 


- 1 


0.79 ± 0.06 


n o/i 
U.y4 


1 

± 


U. id 


Fe 11(2382.77 A)^ 




0.9275 


3.01x10' 




1.01 ± 0.08 


1.21 


± 


0.17 


Al 111(1862.79 A) 


4614.50 


1.4772 


2.79x10" 


1 

-J. 


0.69 ± 0.08 








Zn II('2026 14 A) 

fill XXI \J \J . X JI ± \. ) 


5018 95 


1.4771 


4 89x10" 


-1 


71 ± 07 


1.03 


± 


0.13 


+Mg 1(2026.48 A) b 






1.15x10" 


-1 










Cr 11(2056.25 A) 


5093.54 


1.4771 


1.05x10" 


-1 


0.3 ± 0.04 


0.45 


± 


0.08 


Zn 11(2062.66 A) 


5109.62 


1.4772 


2.56x10" 


-1 


0.7 ± 0.04 


0.52 


± 


0.09 


+Cr 11(2062.23 A)"^ 






7.80x10" 


-2 










Fc 11(2382.77 A) 


5137.88 


1.1563 


3.01x10" 


-1 


0.64 ± 0.08 


0.81 


± 


0.11 


Mg 11(2796.35 A) 


5390.80 


0.9278 


6.12x10" 


-1 


1.04 ± 0.10 


1.06 


± 


0.13 


Mg 11(2803.53 A) 


5403.80 


0.9275 


3.05x10' 


-1 


0.71 ± 0.07 


0.91 


± 


0.08 


Mg 1(2852.96 A) 


5493.01 


0.9254 


1.83 




0.67 ± 0.13 


0.63 


± 


0.15 


Fe 11(2600.17 A) 


5605.97 


1.1560 


2.24x10" 


-1 


1.32 ± 0.07 


1.21 


± 


0.11 


Fe 11(2344.21 A) 


5807.30 


1.4773 


1.10x10" 


-1 


1.97 ± 0.06 


1.86 


± 


0.06 


Fe 11(2374.46 A) 


5881.78 


1.4771 


3.26x10" 


-2 


1.67 ± 0.08 


1.80 


± 


0.22 


Fe 11(2382.77 A) 


5902.60 


1.4772 


3.01x10" 


-1 


2.37 ± 0.09 


2.10 


± 


0.12 


Mg 11(2796.35 A) 


6030.31 


1.1565 


6.12x10" 


-1 


2.07 ± 0.13 


2.08 


± 


0.16 
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Table 1. LRIS Line Identifications (Continued) 



Line(Ava^(A)) 


Ahelio(A) 


z 


fij 








Mg 11(2803.53 A) 


6044.97 


1.1562 


3.05x10" 


-1 


1.77 ± 0.14 


1.79 ± 0.15 


Mn 11(2576.88 A) 


6383.63 


1.4773 


3.51x10" 


-1 


0.62 ± 0.07 


0.50 ± 0.10 


Fe 11(2586.65 A) 


6407.39 


1.4771 


6.84x10" 


-2 


1.53 ± 0.08 


1.5 ± 0.15 


Fe 11(2600.17 A) 


6441.14 


1.4772 


2.24x10" 


-1 


2.48 ± 0.08 


2.2 ± 0.13 


Mg 11(2796.35 A) 


6926.28 


1.4769 


6.12x10" 


-1 


3.06 ± 0.04 


3.00 ± 0.11 


Mg 11(2803.53 A) 


6944.62 


1.4771 


3.05x10" 


-1 


2.88 ± 0.04 


2.90 ± 0.09 


Mg 1(2852.96 A) 


7067.07 


1.4771 


1.83 




1.22 ± 0.05 


1.31 ± 0.09 



^Alternative identification to previous entry. 

''Doublet or blend between lines. 

■^Rest equivalent width on 2001 Feb 22.66. 

^Rest equivalent width on 2001 Feb 23.66. 
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Tablc 2. ESI Line Identifications 



Line(Avac(A)) 


Ahelio(A) 


z 


fij 






Al 11(1670.79 A) 


4138.41 


1.47692 


1.83 




1.87 ± 0.43 


Si 11(1808.01 A) 


4478.18 


1.47686 


2.18x10" 


Q 

-o 


0.76 ± 0.1 


Al 111(1854.72 A) 


4592.39 


1.47606 


5.60x10" 


- 1 


0.41 ± 0.20 


Fe 11(2382.77 A)^ 


4592.39 


0.92733 


3.01x10" 


1 

- i 


0.53 ± 0.26 


Al 111(1854.72 A) 


4594.04 


1.47695 


5.60x10" 


1 

- i 


0.80 ± 0.18 


Fe 11(2382.77 A)*^ 


4594.04 


0.92803 


3.01x10" 


1 


1.03 ± 0.23 


Al III('1862 79 A) 


4614.18 


1.47703 


2 79x10' 


-1 


49 ± 07 


Fe 11(2600.17 A) 


5011.47 


0.92736 


2.24x10" 


-1 


0.41 ± 0.11 


Zn 11(2026.14 A) 


5018.98 


1.47711 


4.89x10" 


-1 


0.91 ± 0.12 


+Mg 1(2026.48 A)^ 






1.15x10" 


-1 




Cr 11(2056.25 A) 


5093.10 


1.47689 


1.05x10" 


-1 


0.44 ± 0.08 


Zn 11(2062.66 A) 


5108.62 


1.47671 


2.56x10" 


-1 


0.69 ± 0.16 


+Cr 11(2062.23 A)'' 






7.80x10" 


-2 




Cr 11(2066.16 A) 


5117.54 


1.47684 


5.15x10" 


-2 


0.27 ± 0.12 


Fe 11(2382.77 A) 


5138.48 


1.15652 


3.01x10" 


-1 


1.12 ± 0.11 


Mg 11(2796.35 A) 


5389.83 


0.92745 


6.12x10" 


-1 


0.96 ± 0.14 


Mg 11(2803.53 A) 


5404.04 


0.92758 


3.05x10" 


-1 


1.08 ± 0.16 


Mg 1(2852.96 A) 


5498.99 


0.92747 


1.83 




0.51 ± 0.21 


Fe 11(2260.78 A) 


5597.47 


1.47590 


2.44x10" 


-3 


0.10 ± 0.08 


Fe 11(2260.78 A) 


5599.68 


1.47688 


2.44x10" 


-3 


0.42 ± 0.16 


Fe 11(2600.17 A) 


5606.38 


1.15616 


2.24x10" 


-1 


1.15 ± 0.16 
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Table 2. ESI Line Identification (Continued) 



Line(Ava,;(A)) 


Ahelio(A) 


z 


fij 






Fe 11(2344.21 A) 


5804.03 


1.47590 


1.10x10' 


-1 


0.37 ± 0.11 


Fe 11(2344.21 A) 


5806.33 


1.47688 


1.10x10" 


_i 


1.51 ± 0.23 


Fe 11(2374.46 A) 


5878.93 


1.47590 


3.26x10" 


_2 


0.30 ± 0.06 


Fe 11(2374.46 A) 


5881.25 


1.47688 


3.26x10" 


_2 


1.24 ± 0.10 


Fe 11(2382.77 A) 


5899.50 


1.47590 


3.01x 10" 


-1 


0.56 ± 0.09 


Fe 11(2382.77 A) 


5901.84 


1 A ''7 o o 

1.47688 


3.01x 10" 


-1 


2.21 ± 0.16 


Mg 11(2796.35 A) 


6028.98 


1.15602 


6.12x 10" 


-I 


2.49 ± 0.08 


Mg 11(2803.53 A) 


6044.44 


1.15601 


3.05x 10" 


-I 


2.14 ± 0.60 


Mg 1(2852.96 A) 


6152.90 


1.15667 


1.83 




0.54 ± 0.31 


Mn 11(2576.88 A) 


6382.58 


1.47686 


3.51x10" 


-1 


0.58 ± 0.12 


Fe 11(2586.65 A) 


6404.29 


1.47590 


6.84x10" 


-2 


0.64 ± 0.13 


Fe 11(2586.65 A) 


6406.82 


1.47688 


6.84x10" 


-2 


1.72 ± 0.20 


Mn 11(2594.50 A) 


6426.18 


1.47685 


2.71x10" 


-1 


0.68 ± 0.12 


Fc 11(2600.17 A) 


6437.76 


1.47590 


2.24x10" 


-1 


0.49 ± 0.05 


Fe 11(2600.17 A) 


6440.31 


1.47688 


2.24x10" 


-1 


2.12 ± 0.19 


Mn 11(2606.46 A) 


6455.96 


1.47691 


1.93x10" 


-1 


0.56 ± 0.12 


Mg 11(2796.35 A) 


6923.66 


1.47596 


6.12x10" 


-1 


1.63 ± 0.24 


Mg 11(2796.35 A) 


6926.76 


1.47707 


6.12x10" 


-1 


1.51 ± 0.19 


Mg 11(2803.53 A) 


6940.97 


1.47580 


3.05x10" 


-1 


1.39 ± 0.15 


Mg 11(2803.53 A) 


6944.06 


1.47690 


3.05x10" 


-1 


1.61 ± 0.11 


Mg 1(2852.96 A) 


7063.74 


1.47593 


1.83 




0.42 ± 0.24 


Mg 1(2852.96 A) 


7066.33 


1.47684 


1.83 




1.21 ± 0.22 



'^Alternative identification to previous entry. 
''Doublet or blend between lines. 
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Table 3. Derived Column Densities for the zu, System 



lon(j) 


Avac(A) 


fij 


b (km s 1) 


log(Ar^) 


log(7VHi)" 


log(iVHl)corr'' 


[Zn/JY 


Mn+ 


2576.88 


0.3508 


22.8 


14.00 ± 0.44 


20.47 ± 0.44 


21.92 ± 0.44 


0.50 




2594.50 


0.2710 


35.3 


13.95 ± 0.17 


20.42 ± 0.17 


21.87 ± 0.18 






2606.46 


0.1927 


29.6 


14.00 ± 0.20 


20.47 ± 0.20 


21.92 ± 0.21 




Mg+ 


2803.53 


3.05x10-1 




> 13.88'! 


> 18.30 


> 19.85 


< 2.65 


Fe+ 


2382.77 


3.01x10-1 




> 14.16^ 


> 18.65 


> 20.92 


< 2.3 



^Assuming solar abundances and no depletion. 

^Assuming depletion corrections from Savage & Sembach (1996). 
^Average log metal abundance ratio relative to Solar value. 
'^Using Equation (3). 
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Tabic 4. Cohimn Densities for Blended Lines 



lon(j) 


log(iVj) 


log(iVHi)^ 


log(iVHl)corr'' 


[Zn/J]^ 


Zn+ 


13.60 ± 0.07 


21.04 ± 0.07 


21.71 ± 0.07 




Cr+ 


14.02 ± 0.1 


20.37 ± 0.1 


22.65 ± 0.1 


0.61 



^Assuming solar abundances and no depletion. 

^Assuming depletion corrections from Savage & Sembach 
(1996). 

'^Log metal abundance ratio relative to Solar value. 
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Table 5. Photoionization Model Parameters 



7VHi(cm-2) 


nHi(cm 3) 


rmax(pc) 


Model 


]^g20.44 


0.089 


1000.0 


I 


]^q20.44 


0.89 


100.0 


II 


-|^g20.44 


89 


1.0 


III 


2Q20.44 


2.75 X 10^ 


6.48 X 10-2 


IV 


"1^022.00 


3.2 


1000.0 


V 


|g22.00 


2.0 X 10^ 


0.19 


VI 


2Q22.00 


10^ 


6.48 X 10-2 


VII 


-|^g22.00 


10^ 


3.56 X 10-2 


VIII 
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Fig. 1.— LRIS spectra of GRB 010222 on 2001 Feb. 22.66 and Feb 23.66 corrected for 
Galactic foreground extinction and smoothed with a boxcar corresponding to the instru- 
mental resolution. Three absorption systems are labelled: zi — 1.47688 {dashed lines), 
zi — 1.15628 {solid lines), and — 0.92747 {dotted lines). Letters in parentheses refer to 
the notes in Table 1. 
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Fig. 2.— Continuum-normalized ESI spectrum of GRB 010222 on 2001 Feb. 23.61. The 
spectrum is smoothed with a boxcar corresponding to the instrumental resolution. Three 
absorption systems are labelled: zi — 1.47688 {dashed lines), zi — 1.15628 {solid lines), and 
zs — 0.92747 {dotted lines). Atmospheric bands are also indicated. Letters in parentheses 
refer to the notes in Table 2. 
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Fig. 3. — Optical continuum points with la standard deviation uncertainties, along with the 
best fitted line that minimizes for 2001 Feb 22.66 and 2001 Feb 23.66. 
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Fig. 4. — Multiwavelenth X-ray and optical continuum corrected for dust using the extinction 
curve of Calzetti et al. (2000) and color excess E{B - V) ^ 0.1 for 2001 Feb 22.66. The 
resulting optical spectrum has a slope P —0.43. The X-ray data is shown schematically 
as a dashed line, while the best fit through the optical data is plotted as a solid line. 
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Fig. 5. — Metallic absorption-line profiles in the ESI spectrum plotted in velocity space, 
showing evidence of substructure in the most distant absorber (2:1). As zero velocity we use 
the systemic redshift, z — 1.47688. The dashed line indicates the average offset of = —119 
km s~^ from the systemic velocity. 
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Fig. 6.— Modelled time evolution of the Mg I (2852.96 A) and Mg II (2796.35 A) equivalent 
widths for different sets of initial conditions (Table 5), constrained by (A^hi = lO^^'^"^ cm~^). 
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Fig. 7.— Modelled time evolution of the Mg I (2852.96 A) and Mg 1(2796.35 A) equivalent 
widths for different sets of initial conditions (Table 5), constrained by (A^hi = 10^^ cm~^). 



